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Abstract 

An elastic string embedded between rigid walls is simulated by means of the density-matrix 
renormalization group. The string collides against the walls owing to the quantum- mechanical 
zero-point fluctuations. Such "quantum entropic" interaction has come under thorough theoretical 
investigation in the context of the stripe phase observed experimentally in doped cuprates. We 
found that the excitation gap opens in the form of exponential singularity AE ~ exp(— Ad a ) {d: 
wall spacing) with the exponent a = 0.6(3), which is substantially smaller than the meanfield 
value a = 2. That is, the excitation gap is much larger than that anticipated from meanfield, 
suggesting that the string is subjected to robust pinning potential due to the quantum collisions. 
This feature supports Zaanen's "order out of disorder" mechanism which would be responsible to 
the stabilization of the stripe phase. 

PACS numbers: 74.72.-h High-Tc compounds, 05. 10. -a Computational methods in statistical physics and 
nonlinear dynamics, 46.70.Hg Membranes, rods and strings, 05.70.Jk Critical point phenomena. 
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I. INTRODUCTION 



Recently, Zaanen brought up a problem of "quantum string," which is a linelike object 
subjected to line tension, and it wanders owing to quantum zero-point fluctuations^!. Cen- 
tral concern is to estimate the interaction among strings as it wanders quantum-mechanically 
and undergoes entropy-reducing collisions with adjacent neighbors. Statistical mechanics of 
quantum-string gas would be responsible to the low-energy physics of the stripe phase ob- 
served experimentally in doped cupratesSBi. In particular, one is motivated to gain insights 
how the stripe pattern acquires stability. Actually, a good deal of theoretical analyses had 
predicted tendency toward stripe-pattern formation!^!'!!. However, ab initio simulations on 
the t-J model still remain controversial about that issue0'00l!0. 

In the series of papersB'i, the authors pointed out that the "quantum entropic" interaction 
lies out of conventional picture, and it would rather give rise to the solidification of the 
gas of strings; namely, the stripe phase is stabilized by mutual collisions. (In a different 
context, the entropic interaction was explored in Refs. |15|,16|. ) Their analysisSS is based 
on the Helfrich approximation^, which has been very successful in the course of studies of 
stacked membranes under thermal undulations. Based on this approximation, they reveled 
a significance of /ong^wavelength fluctuations. Note that in the conventional picture, on 
the contrary, the string is "disordered" as in the Einstein-like view of crystal, and only 
short-wavelength fluctuations contribute to collisions. 

To be specific, the key relation in their theory is the following, 

with collision-induced energy cost / (per unit volume) and elasticity modulus (with respect 
to the compression of string intervals) B = d 2 d 2 f jdd 2 \ see the Hamiltonian @ as well. The 
logarithmic term signals the significance of long-wavelength fluctuations. The relation ([!]) 
yields remarkable consequences. For instance, the elastic modulus is given by the stretched 
exponential form B ~ exp It is noteworthy that the elasticity modulus B is much 

larger than that of meanfield B ~ exp(— Ad 2 ). Hence, the result indicates that contrary to 
our naive expectations, short-wavelength modes are suppressed owing to collisions, and the 
collisions rather contribute to the solidification of the string gas; namely, the stripe phase is 
stabilized by the "order out of disorder" mechanismBi. 
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In our preceding papeiO, we have verified the above relation (HJ), and demonstrated 
that the elastic modulus is actually governed by the stretched exponential form B ~ 
exp(— Ad ' 808 ^ 1 '). We have carried out the first-principle simulation by means of the density- 
matrix renormalization grouplil00: We put a quantum string between rigid walls with 
spacing d, and measured its repelling interaction. (This technique has been utilized in the 
course of studies of fluctuation pressure of stacked membranesH'fiH. ) The Hamiltonian is 
given by, 

n = il(^- + vixi)) + £ §(*« - (2) 

i=l \ / 1=1 

Here, denotes the operator of transverse displacement of a particle at zth site, and pi is 
its conjugate momentum. They satisfy the canonical commutation relations [a^Pj] = ih5y, 
= and \pi,Pj] = 0. V(x) is the rigid-wall potential with spacing d; 

\ for < x < d 
V(x) = { ~ ~ . (3) 

I oo otherwise 

E denotes the line tension which puts particles into line. Classical version of this Hamiltonian 
has been used as a model for line dislocations and steps on (vicinal) surfaces§lS0. Note 
that for sufficiently large S, one can take continuum limit, with which one arrives at field- 
theoretical version of quantum string^. 

In the present paper, we are concerned in the excitation (mass) gap due to the collisions. 
We postulate the following exponential singularity, 

AE ~ exp(-AT), (4) 

for sufficiently large d. In Ref. 0, the gap (crossover temperature) T is calculated in the 
form To oc VB; namely, one should obtain a = 2/3. This result cannot be understood in 
terms of the meanfield picture yielding much smaller mass gap a = 2. (We will outline this 
picture in the next section.) The purpose of this paper is to judge the validity of those 
scenarios by performing first-principle simulations. Our result is a ~ 0.6(3) for S = 0.5-4. 

The rest of this paper is organized as follows. In the next section, we argue physical 
implications of the singularity exponent a. We introduce the viewpoint of critical phenomena 
in order to interpret the gap formula (HI). The meanfield argument is also explicated. In 



Section III, we perform numerical simulations, and show evidences of the breakdown of 



meanfield. In the last section, we give summary and discussions. 



II. INTERPRETATION OF THE SINGULARITY EXPONENT a — A 

STATISTICAL-MECHANICAL OVERVIEW 



In Introduction, we have overviewed the quantum string in the context of the stripe phase. 
In this section, we will survey different aspects of the quantum string in terms of statistical 
mechanics. This viewpoint is utilized in the simulation-data analyses in the succeeding 
section [Enj. 

In the path-integral picture, a quantum string spans a "world sheet" as time evolves^. 
Therefore, a quantum string is equivalent to the random surface under thermal undulations. 
The random surface is in the critical phase of Kosterlitz and Thouless. (In our model @, 
however, there appears no "flat phase," because the "height" variables {x{\ are continuous.) 
The presence of rigid walls is supposed to destroy the KT criticality. Therefore, the collision- 
induced mass gap (^) reflects the universality class of the phase transition at 1/d — > 0. In 
that sense, our aim is to determine the universality class. 

The phase transition has not been studied very extensively so far. The rigid-wall potential 
is, by nature, non-analytic, and hence, it is quite cumbersome to carry out perturbative 
analyses such as loop expansions. Moreover, the confinement due to walls destroys the 
applicability of the duality transformation, with which one could establish an elucidating 
equivalence between the random surface and the two-dimensional XY model. 

Even in computer simulation, the verification of the gap formula (^) is rather troublesome. 
Note that the gap opens very slowly. That is, the correlation length is kept exceedingly large 
in the vicinity of the critical point. The correlation length would exceed the system sizes. 
Note that the conventional criticality exhibits the power law singularity, 

AE ~ ~ (l/dy, (5) 

which would be much easier to cope with. We will overcome the difficulty through resorting 
to the density-matrix renormalization grou ; with this technique, we treat system sizes 

up to L = 42. 

Now, we are in a position to discuss the meaning of the singularity exponent a. The 
correlation length is proportional to the inverse of the excitation (mass) gap; namely, £ ~ 
exp(Acf"). £ sets the characteristic length scale of the collision intervals. Therefore, in 
the length scale I = £, the random surface fluctuates freely. And so, the mean fluctuation 
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deviation is calculated 



^-(log/) 1 / 2 -^ 2 . (6) 
From this relation, we see that a reflects an amount of fluctuation amplitudes. In the 



meantime, the meanfield picture insists J (xj) oc d. (More precisely, an artificial mass term 
is included so as to enforce the condition; see Refs. |16|,^2[) Hence, one should obtain, 



a = 2, (7) 

in the meanfield picture. 

On the contrary, as is mentioned in Introduction, recent theoryBi predicts much smaller 



singularity exponent a = 2/3, from which one obtains a counterintuitive result y(x 2 ) ~ 
d 1 / 3 <C d. The result indicates that the string is straightened macroscopically by collisions: 
This is quite reminiscent of the aforementioned claim that the infrared modes are still active 
withstanding the collisions^!. Similar, but slightly modest, exponent a = 1 was reported 
(in a different context) with analytical calculations^!^. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we perform first-principle simulations so as to estimate the singularity 
exponent a. First, we will explicate the simulation algorithm. 

A. Details of the density-matrix renormalization group 

The density-matrix renormalization group has been utilized for the purpose of treating 
very large system size The technique is based on an elaborate reduction of the Hilbert- 
space bases, and it has been applied to various spin and electron systems very successfully^. 
Meanwhile, the technique had become applicable to bosonic systems such as phonons and 
oscillat or sSiHiHH. (Note that the full exact diagonalization does not apply, because even a 
single oscillator spans infinite-dimensional Fock space. The huge dimensionality overwhelms 
the computer memory space.) Here, we employ the density- matrix-renormalization tech- 
nique in order to diagonalize the quantum string (Q). Full account of technical details is 
presented in our preceding papei@. In this paper, we will outline the simulation algorithm 
with an emphasis on the modification to calculate the excitation gap precisely. 
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To begin with, we need to set up local (on site) Hilbert-space bases: Provided that the 
line tension £ is turned off, oscillators are independent, and each of them reduces to the 
text-book problem of "particle in a box." Hence, the full set of eigensystems is calculated 
easily. In this way, we set up the intra-site Hilbert space with use of the low-lying M — 12 
states, and discarded the others. 

Secondly, provided that those intra-site bases are at hand, we are able to apply the 
density-matrix renormalization group. The total system consists of block, site, site and 
block. The left-half part (block and site) is renormalized into a new block with reduced 
dimensionality m; we set, at most, m = 50. More precisely, those m bases are chosen from 
the eigenvectors of the density matrix for the left-half part, 

p = Tr R |*><*|, (8) 

with the largest m eigenvalues. Tr R denotes the partial trace with respect to the right-half 
part, and = |0) + |1) with the ground state |0) and the first excitation state |1) of the 
total system. Note that in our previous simulationS, we simply set \^f) = |0), because we 
were just concerned in the ground state properties. 

To summarize, we had carried out truncations of bases through two steps: One is the 
truncation of the intra-site bases M, and the other is the truncation of the "block" states m 
through the density-matrix renormalization. Each of these procedures is monitored carefully, 
and the relative error of low-lying energy levels is kept within 10~ 7 ; see Ref. [TBI for details. 



We have repeated renormalizations 20-times. The system size reaches L = 42. 



B. Excitation energy gap 

We have plotted the energy gap against the inverse of wall spacing 1/d for several values 
of line tension E; see Fig. |IJ From the plot, we see that the gap opens extremely slowly; 
namely, in the range 1/d < 0.3, the mass gap AE is maintained to be very small. This is 
a typical signature of the exponential singularity (|j) rather than the power law (|5|). For 
1/d < 0.2, the simulation does not continue, because the numerical diagonalization fails in 
resolving the nearly degenerate low-lying levels. 

In the inset, we have presented the logarithmic plot. The data exhibit convex curves. 
Hence, we see that the singularity exponent a would not exceed a = 1. As is mentioned in 
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Section ||, the meanfield argument predicts a = 2. Therefore, we found that the collision- 
induced gap cannot be understood in terms of the meanfield picture. 



C. j3- function analyses 

The criticality is best analyzed by the /^-function. The /3-function describes the flow of 
a certain controllable parameter (in our case, 5 = 1/d) with respect to the infinitesimal 
rescaling of the unit of length; namely, 

with correlation length £. Because the mass gap is the inverse of the correlation length, we 
obtain, 

« 1 /* = -(d^o b,AB )" 1 ' (10) 

with excitation gap AE. 

As is emphasized in Section [TT], the gap formula ([|) is deeply concerned with the criticality 
at 5 — 1/d — *■ 0. The /5-function reflects the universality class of the phase transition. For 
instance, as for the exponential singularity such as Eq. it behaves like, 

p(i/d) ~ \/{Ao){\/dy + \ (ii) 

On the other hand, the conventional second-order transition (|5]) is characterized by the 
behavior, 

0(l/d) = (l/d)/u. (12) 

In this way, we can read off the critical exponents such as a and v. 

As is presented in the previous subsection, the excitation-gap data AE are readily avail- 
able. Remaining task is to perform the numerical derivatives of AE'. We adopted the 
"Richardson's deferred approach to the limit" algorithm in the text book0. In this algo- 
rithm, one takes an extrapolation after calculating various finite-difference differentiations. 
We monitored the relative error, and checked that the error is kept within 10~ 8 . From these 
preparations, we can calculate the /3-function from our simulation data. 

We plotted the /9-function (3(1/ d) in Fig. || for S = 2. We see that there appear two 
regimes: As is indicated in the plot, in the region 1/d ^ 0.3, the /^-function is best fitted by a 
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power law (l/cf) L6 . Hence, the universality would belong to the exponential singularity ((TTJ) 



with exponent a = 0.6. Note that such stretched-exponential behavior is suggested in the 
previous subsection. On the contrary, the /3-function falls in the simple behavior l/v{l/d) 



box." Hence, the excitation gap is simply given by the power law (j|) with the index v = 2. 
To summarize, we found that for large 1/d, the physics changes so that we must look into 
the behavior in the vicinity of 1/d = 0. 

In order to estimate a, we have presented the log- log plot in Fig. |3| for various S. From 
the slopes, we obtain an estimate o = 0.6(3). This estimate is far less than the mean-field 
value a = 2. It clearly supports Zaanen's claim a = 2/3011 that the collision-induced mass 
gap is governed by the stretched-exponential singularity. That is, the collision-induced mass 
gap is far larger than that anticipated from meanfield, and hence, the string is subjected to 
robust pinning potential due to quantum collisions; see Section |I[ 

Meanwhile, we had also calculated the Roomany-Wyld approximant of the /3-functionil, 



where / and /' denote a pair of system sizes. This approximant, in general, exhibits smaller 
finite-size collections. In fact, we found that it converges rapidly to the thermodynamic 
limit. However, the resultant data are identical to those with the aforementioned formula 
(|l0l) for sufficiently large L. Because we take advantage of very large system sizes owing to 
the density-matrix renormalization group, the approximant formula (|i~3|) is not particularly 
necessary. 

D. Mean- fluctuation deviation 

So far, we have extracted the singularity exponent directly from the excitation-gap data. 
However, as is mentioned in Section |J, the mean deviation of the string undulations contains 
an information of the exponent; 




1 + ]n(AE l /AE v )/ln(l/l') 

{KeJKeJJae^Ke^T 2 



(13) 




(14) 
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In Fig. |], we plotted the mean fluctuation width A for various d and £ = 0.5-4. From 
the plot, we see that there appear two distinctive regimes as in Fig. |2|. For small d, A is 
proportional to d, as is anticipated from a naive argument; see Section [TT[ However, for 
large d, A starts to bend down, indicating that the string undulation becomes suppressed. 
Hence, for large d, the conventional picture does not apply. From the slopes, we read off 
the exponent a/2 « 0.5. (This result supports our estimate a = 0.6(3) rather than the 
meanfield value a = 2.) Hence, we see A/d <C 1 as for d — ► oo. As a consequence, we 
found fairly definitely that the string is straightened macroscopically due to the quantum 
collisions. This feature is reminiscent of the infrared divergences encountered in the theory 
by Zaaneni. 



IV. SUMMARY AND DISCUSSIONS 



We have investigated the quantum entropic (collision induced) interaction of an elastic 
string confined between walls (0). We performed the first-principle simulation by means 
of the density- matrix renormalization group. Our aim is to estimate the singularity expo- 
nent er in the excitation-gap formula (f|). The exponent contains a number of significant 
informations which are overviewed in Section |H|. 

In order to estimate a, we utilized the /^-function (|T0|) , which is readily accessible from our 
simulation data. From the asymptotic form of large d, we obtained the estimate a = 0.6(3) 
for £ = 0.5-4; see Fig. §. Our estimate supports the recent proposafli that AE is described 
by the stretched exponential form with a = 2/3, whereas it might conflict with the meanfield 
value <j = 2. It would be noteworthy that the elastic modulus also obeys the stretched 
exponential, as is shown in our previous studyEl. Hence, it is established fairly definitely 
that the excitation gap is also described by the stretched exponential singularity. Namely, 
the excitation gap is much larger than that anticipated from meanfield. Moreover, we found 
that the mean fluctuation width A is far less than the wall spacing d; namely A ~ d a5 . This 
again supports our conclusion rather than the meanfield; see Eq. (|5|). 

From those observations, we are led to the conclusion that the string is confined by 
extremely robust effective pinning potential due to the quantum collisions; subsequently, 
large mass gap opens and A gets bounded. Those features are also quite consistent with 
the findings of the series of worksEH, validating their "order out of disorder" mechanism 
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FIG. 1: Excitation gap AE is plotted for 1/d (d: wall spacing). The excitation gap AE opens 
extremely slowly: This feature is characteristic of the universality class of exponential singularity 
(Q) rather than the power law (||). The inset shows the logarithmic plot. The curves bend convexly. 
Therefore, the singularity exponent would not exceed a = 1. 

responsible to the stabilization of the stripe phase. 

Provided that plural strings are concerned, does the value a get affected? As a matter of 
fact, the exponent a is not necessarily universal, but it does exhibit continuous variationB^I 
for a certain class of extended models. More specifically, in Refs. the authors claim 

that for extended models with high central charges (that is, many random surfaces) exhibit 
various types of criticalities. Therefore, It is likely that a would acquire corrections for 
multi-string case. Extended simulation directed to this issue would be remained for future 
study. 
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FIG. 2: /3-function ( |T0| ) is plotted for S = 2. There appear two regimes with respective asymptotic 
forms; see text. The small-l/cf behavior indicates that the transition belongs to the universality 
class of exponential singularity with a ~ 0.6; see Eqs. (11) and (|l^). 
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